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We report first principles density functional calculations for 5,6-diliydroxyindole-2-carboxylic acid 
(DHICA) and several reduced forms. DHICA and 5,6-diliydroxyindole (DHI) are believed to be the 
basic building blocks of the eumelanins. Our results show that carboxylation has a significant effect 
on the physical properties of the molecules. In particular, the relative stabilities and the HOMO- 
LUMO gaps (calculated with the ASCF method) of the various redox forms are strongly affected. 
We predict that, in contrast to DHI, the density of unpaired electrons, and hence the ESR signal, 
in DHICA is negligibly small. 

PACS numbers: 87.15.-v, 82.35.Cd, 87.64.Je 



I. INTRODUCTION 

The melanins are an important class of pigmentary 
macromolecule found throughout nature from fungi to 
man Two types of melanin are found in the skin, 
hair, eyes and nervous systems of humans: pheome- 
lanin and eumelanin. Pheomelanin is associated with 
ginger hair and is a cysteinyl-dopa derivative. Eume- 
lanin, which we consider here, is found in much higher 
levels in black hair and is believed to be a macromolecule 
of 5,6-dihydroxyindole-2-carboxylic acid (DHICA) and 
5, 6-dihydroxy indole (DHI) 1, 2]. Perhaps the most 
surprising property of melanins is that, whereas most 
biomolecules show only well defined absorption peaks 
around 280 nm (4.5 eV) and little absorption at lower 
energies 4], the melanins show broad band monotonic 
absorption in the range 1.5 to 5 eV Q. Melanins are 
also efficient free radical scavengers and antioxidants Q| . 
Thus melanins are ideally suited for their best understood 
biological function as the primary photoprotectant in our 
skin and eyes. Conversely, both pheomelanin and eu- 
melanin are implicated in the development of melanoma 
skin cancer 5[. Melanins also display semiconductor like 
behaviour 0^3, that suggests possible uses in applica- 
tions such as bio-sensors and bio-mimetic photovoltaics 
1^, However, it is not clear that the conductivity re- 
ported in any of these experimental studies is electronic 
in nature |il|Tl|. 

Surprisingly, given the obvious importance of the 
melanins and several decades of work, little is known 
about the general structure-property-function relation- 
ships that control their behaviour. This is, at least in 
part, due to the difficulties in studying the melanins; they 
are chemically and photochemically stable, and are vir- 
tually insoluble in most common solvents. For example, 
major questions still remain concerning their basic struc- 
tural unit 1^. Two opposing schools of thought exist: 
i) that eumelanins are composed of highly cross-linked 
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extended hetero-polymers based upon the Raper-Mason 
scheme Ij, and ii) that eumelanins are actually composed 
of much smaller oligomers condensed into 4 or 5 oligomer 
nano-aggregates l^J. A clear idea of the basic structural 
unit is critical to developing a consistent model for con- 
densed phase charge transport in such disordered organic 
systems. It is also important in the context of "molecu- 
larly engineering" melanins to have the ability to create 
or enhance functionality in high technology applications 

Recently the role of disorder in determining the phys- 
ical properties of the melanins has been emphasised 
I14L Il5| . In particular it has been proposed that 
the broad band monotonic optical absorption of eume- 
lanin arises from the existence of a large number of dif- 
ferent monomers and hence a wide variety of macro- 
molecules. The idea is, basically, that, although each 
of these molecules has a distinct optical absorption that 
is not atypical for bio-macromolecule, summing over a 
large range of absorption frequencies in a biological en- 
vironment (in which both dipole fluctuations jl^ and 
thermal effects broaden the spectrum) gives rise to the 
observed broad band monotonic absorption spectrum. 

Despite the widespread agreement that natural eume- 
lanin contains both DHI and DHICA, previous quantum 
chemical studies 14, il [ll [ll S 111 [23 ] have 
only considered DHL This is in part due to the greater 
chemical complexity of DHICA, but also in part due to 
assumption that the carboxylation of DHI does not play 
a significantly role in determining the properties eume- 
lanin. However, given the potential importance of the 
role of disorder in the physical properties of eumelanin 
fidi Q, a careful study of the similarities and dif- 
ferences between DHI and DHICA is clearly required. 
Our study is also motivated by recent work on synthetic 
DHICA and the need for high quality first prin- 
ciples calculations to compare with future experiments 
on synthetic DHICA. In this paper we present, to our 
knowledge, the first quantum chemical study of DHICA. 
We use density functional theory (DFT) to calculate the 
physical and electronic structure of DHICA and several 
reduced forms, each of which is drawn schematically in 




FIG. 1: Schematic representation of 5,6-dihydroxyindole-2- 
carboxylic acid (DHICA) 1 and the reduced forms which we 
also consider here. The numbers correspond to those used in 
table H] Note that 2a and 2b are radicals in the neutral state. 

figure n Our results show that there are significant 
differences between the physical properties of DHI and 
DHICA. Specifically we show that the relative stability 
of the reduced forms is significantly altered by carboxy- 
lation as is the HOMO-LUMO gap [the energy difference 
between the highest occupied molecular orbital (HOMO) 
and the lowest unoccupied molecular orbital (LUMO)]. 

II. CALCULATION DETAILS 

The chemical and electronic structures were found 
from first principles DFT calculations. We have per- 
formed our calculations using the Naval Research Lab- 
oratory Molecular Orbital Library (NRLMOL) [H 
miH, S 113, NRLMOL performs massively par- 
allel electronic structure calculation using gaussian or- 
bital methods. Throughout we have used the Perdew, 
Burke and Ernzerhof (PBE) 32] exchange correlation 
functional, which is a generalised gradient approximation 
(GGA) containing no parameters. For each molecule we 
have fully relaxed the geometry with no symmetry con- 
straints. 



TABLE I: The calculated bond lengths (in A) and bond an- 
gles (in degrees) of DHICA. The atom numbers correspond 
to those shown in figure Q Note that there is no experi- 
mental data or previous theoretical results to compare with. 
We therefore contrast our results with previously published 
results for DHI 15]. 





DHICA 


DHI 


N1-C2 


L387 


1.390 


C2-C3 


L388 


1.365 


C3-C9 


L421 


1.443 


C9-C4 


L414 


1.366 


C4-C5 


L380 


1.458 


C5-C6 


1.428 


1.586 


C5-O10 


1.384 


1.226 


C6-011 


1.365 


1.226 


C6-C7 


1.387 


1.457 


C7-C8 


1.403 


1.361 


C8-C9 


1.432 


1.483 


C8-N1 


1.372 


1.381 


C2-C12 


1.453 


_ 


C12-013 


1.368 


_ 


C12-014 


1.225 


_ 


N1-C2-C3 


109.1 


110.8 


C2-C3-C9 


107.0 


107.2 


C3-C9-C4 


133.5 


132.7 


C3-C9-C8 


107.2 


106.2 


C9-C4-C5 


118.3 


119.7 


C4-C5-C6 


121.7 


118.3 


C4-C5-O10 


124.8 


122.7 


O10-C5-C6 


113.6 


118.8 


C5-C6-011 


118.9 


119.1 


C5-C6-C7 


121.2 


117.9 


011-C6-C7 


119.8 


122.8 


C6-C7-C8 


117.2 


118.3 


C7-C8-C9 


122.3 


124.6 


C8-C9-C4 


119.3 


121.0 


N1-C8-C9 


107.3 


106.0 


C7-C8-N1 


130.5 


129.3 


C8-N1-C2 


109.4 


109.6 


N1-C2-C12 


118.9 




C3-C2-C12 


132.0 




C2-C12-013 


112.9 




C2-C12-014 


124.2 




013-C12-014 


112.8 





III. RESULTS AND DISCUSSION 

We report the calculated bond lengths and bond an- 
gles of DHICA in tabled For comparison we also report 
the same quantities for DHI from calculations 15] using 
the same methods (it should be noted that the results for 
DHI are in good agreement with other quantum chemi- 
cal calculations ^l|, [23, 12^ ) • Unsurprisingly there is lit- 
tle difference between the physical structure of DHI and 
DHICA. 

In tabled we report the calculated energies of each of 
the molecules considered here in there neutral, -1 and -2 
charge states. We also report the relative concentration 
at 300 K in table Hill It can be seen that in the neutral 



state DHICA is considerably more stable than any of the 
reduced forms. Whereas in the -1 state 2a and 2b are 
most stable. This result has important implications for 
electron spin resonance (ESR) experiments. ESR exper- 
iments are sensitive to unpaired electrons which, these 
calculations indicate, are not found in DHICA (this is 
not the case for DHI 0, [2^ ) . Given the recent synthe- 
sis of DHICA, ESR experiments on pure DHICA sam- 
ples would provide a direct check of these calculations. 
These results suggest that samples of DHICA have signif- 
icantly less chemical disorder than pure DHI or natural 
eumelanin. We can therefore speculate that the opti- 
cal absorption of macromolecules formed from synthetic 
DHICA may not be a smooth as that of natural eume- 
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TABLE II: The total energy (in eV) of DHICA and the re- 
duced forms shown in figure^in the neutral, -1 and -2 charge 
states. The energies are quoted relative to the energy of neu- 
tral DHICA. In each case we fully relaxed the geometries as 
shown in figure |21 





Neutral 


-1 


-2 


1 





-0.299 


+3.245 


2a 


+ L165 


-L478 


+L887 


2b 


+ L152 


-L613 


+ L813 


3a 


+2.392 


-0.223 


+ L814 


3b 


+2.693 


-0.307 


+ L273 


3c 


+2.583 


-0.321 


+ L276 


TABLE III: 


The relative abundances of DHICA and the re- 


duced forms shown in figure Q in the neutral, -1 and 


-2 charge 


states at 300 K. We have assumed a Boltzmann distribution 


of each charge state separately with an infinite reservoir of 


H2. The relevant energies 


are stated in table UTI 






Neutral 


-1 


-2 


1 


1 


< 10"^^ 


< 10"^^ 


2a 


< 10-^^ 


5.5 X 10"^ 


< 10-^° 


2b 


< 10-^^ 


1 


<io-^ 


3a 


< 10-^° 


< 10"^^ 


< 10-^ 


3b 


< 10-^^ 


< 10-^^ 


1 


3c 


< 10-^^ 


< 10-^^ 


0.88 



lanin. Note that for all of the molecules considered here 
the -1 charge state is more stable than the neutral state 
(as is the case for the related DHI molecules 14, 15] ), this 
is clearly related to free radical scavenging properties of 
eumelanin. 

One of the most important properties of melanin 



TABLE IV: The HOMO-LUMO (highest occupied molecular 
orbital- lowest unoccupied molecular orbital) gap in eV. The 
HOMO-LUMO gap has been calculated from both the ASCF 
method and from a simple interpretation of the Kohn-Sham 
eigenvalues. We have previously shown 10, 15] that for DHI 
and related molecules the value of the HOMO-LUMO gap 
calculated by ASCF method is in good agreement with the 
results of time dependent density functional theory (TDDFT) 
calculations for the same molecules |23|. The figures in brack- 
ets show the results of equivalent calculations for DHI 15]. It 
can be seen that there is a significant change in the HOMO- 
LUMO gap for 1, which is essentially the only stable monomer 
in the neutral state (c.f. tables UTI and IIII|) . 

Simple interpretation of the ASCF 
Kohn-Sham eigenvalues 



1 2.85 (3.48) 3.04 (3.61) 
2a 2.24 (-) 2.67 (-) 

2b 2.37 (-) 2.64 (-) 

3a 0.87 (1.07) 1.96 (2.02) 

3b 0.79 (0.80) 1.10 (1.12) 
3c 0.89 (-) 1.25 (-) 




FIG. 2: (Colour online.) The electron density in the high- 
est occupied molecular orbital (HOMO) (left) and the lowest 
unoccupied molecular orbital LUMO (right) of, from the top 
down, 1, 2a, 3a and 3b. Although the electron density in 
the HOMO of 1 is very similar to the electron density in the 
HOMO of DHI (reported in [l^) there are significant differ- 
ences between the electron densities in the LUMOs of DHICA 
and DHI near atom C2 (where the carboxylic group attaches) . 
This is presumably related to the large difference (^ 20%) in 
HOMO-LUMO gap calculated for the two compounds (c.f. ta- 
ble|I3). Both the HOMO and LUMO electron densities of 3a 
and 3b are very similar to the electron density in the HOMO 
and LUMO of the equivalent reduced forms of DHI (the in- 
dolequinone and the semiquinone respectively, see 15]) This 
is presumably related to the similarity of HOMO-LUMO gap 
calculated for the carboxylated and uncarboxylated forms of 
the two compounds (c.f. table II V) ). 

monomers is the HOMO-LUMO gap. This is indicative 
of the optical spectrum as it represents the fundamen- 
tal optical absorption and Hiickel model calculations 
for DHI show that the HOMO-LUMO gap of a single 
molecule correlates with the semiconducting gap of an in- 
finite homopolymer. The "semiconductor" gap is likely 
to be a critical design parameter for molecularly engi- 
neered forms of eumelanin which may be useful as func- 
tional materials in electronic devices and sensors . 
Thus controlling the HOMO-LUMO gap of monomers 
and thence of macromolecules may provide a route to 
controlling the semiconducting gap of eumelanins [TEj . 

DFT is a theory of the ground state, therefore cal- 
culations represent the energy gap between Kohn-Sham 
eigenvalues and not the true HOMO-LUMO gap of the 
molecules. This is known as the band gap problem 
|3^] . Additionally, it is accepted that the PEE functional 
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can significantly underestimate the HOMO-LUMO gap . 
Therefore we have also employed the ASCF method 
to calculate the HOMO-LUMO gap. We have previously 
shown that equivalent results for DHI reproduce the 
trends found in time dependent DFT calculations [2^. 
In table HSI we compare the HOMO-LUMO gap found 
from a simple interpretation of the Kohn-Sham eigenval- 
ues with those found by the ASCF method. We also re- 
produce the equivalent results for some of the uncarboxy- 
lated forms from From this it can be seen that for 3a 
and 3b carboxylation does not significantly change the 
HOMO-LUMO gap. But, for 1, which our results show 
is the only significant component of DHICA samples {ta,- 
ble IHI|) . carboxylation causes a significant change in the 
HOMO-LUMO gap. This indicates that, because natural 
eumelanin contains both DHI and DHICA, a sample of 
eumelanin contains a large range of HOMO-LUMO gaps, 
even before the effects of the formation of macromolecules 
are considered. 

The changes in the HOMO-LUMO gap in DHI are un- 
derstood in terms of the delocalisation of the wavefunc- 
tion. Therefore in figure [2 we plot the calculated elec- 
tron density of 1, 2a, 3a and 3b respectively. Comparing 
these plots with the equivalent results for the uncarboxy- 
lated forms we see that for 3a and 3b the electron 
densities in both the HOMOs and LUMOs is remarkably 
similar in both the carboxylated (figure Ej) and uncar- 
boxylated forms. This is presumably the reason that the 
calculated HOMO-LUMO gap is essentially independent 
of carboxylation (table II V|) . On the other hand while 
the electron density of the HOMO of 1 is very similar 
to that of its uncarboxylated form, there are significant 
differences between the electron density of the LUMO 
and that of the uncarboxylated form (figure [2|) . This 
is, presumably, associated with the large change in the 
HOMO-LUMO gap caused by carboxylation (c.f. table 
llV|l . 



IV. CONCLUSIONS 

We have carried out first principles density func- 
tional calculations for 5,6-dihydroxyindole-2-carboxylic 



acid (DHICA) and several reduced forms. To the best 
of our knowledge these are the first quantum chemical 
calculations for this system. These molecules and similar 
molecules based on 5,6-dihydroxyindole (DHI) are be- 
lieved to be the basic building blocks of the eumelanins. 
We used the ASCF method to study the HOMO-LUMO 
gap. Comparing our results for DHICA with previously 
published results for DHI we found that that the addi- 
tion of the carboxylic acid group has a significant effect 
on the physical properties of the molecules and in partic- 
ular on the relative stability of the various redox forms 
and the HOMO-LUMO gap. Based on the calculated 
relative stabilities of the various redox forms in their var- 
ious charge states we find predict that there will be an 
extremely low unpaired electron density in samples of 
pure DHICA (this is not the case for DHI HI [13) and 
thus a negligible signal should be seen in ESR experi- 
ments. This indicates that experimental results cannot 
be straightforwardly extrapolated from DHICA to DHI 



or vice versa. 
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